Travelling wave solutions to nonlinear Schrodinger equation with self-steepening and 

self-frequency shift 



O 

(N 



Anshul Saini, 1 Vivek M. Vyas, 1 S. N. Pandcy, 2 T. Soloman Raju, 3 and Prasanta K. Panigrahi 1 

1 Indian Institute of Science Education and Research (USER)- Kolkata, Mohanpur, Nadia-7^1252, India 
2 Motilal Nehru National Institute of Technology, Allahabad - 211004, India 
3 TempleCity Institute of Technology & Engineering, Bhubaneswar 752 057, India 

We investigate exact travelling wave solutions of higher order nonlinear Schrodinger equation in 
the absence of third order dispersion, which exhibit non-trivial self phase modulation. It is shown 
that, the corresponding dynamical equation, governing the evolution of intensity in the femtosecond 
regime, is that of non-linear Schrodinger equation with a source. The exact localized solutions to this 
system can have both super and subluminal propagation belonging to two distinct class. A number 
of these solitons exhibit chirality, thereby showing preferential propagation behavior determined 
by group velocity dispersion. Both localized bright and dark solitons are found in complementary 
velocity and experimental parameter domains, which can exist for anomalous and normal dispersion 
regimes. It is found that, dark solitons in this system propagate with non-zero velocity, unlike their 
counterpart in nanosecond regime. Interestingly, subluminal propagation is observed for solitons 
having a nontrivial Pade type intensity profile. 
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I. INTRODUCTION 

Optical solitons were first predicted by Hasegawa and 
Tappert [lj and later observed by Mollenauer et. al. 
Q. The existence of soliton solutons in non-linear 
Schrodinger equation (NLSE) is related to the complete 
integrability of this dynamical system (3|. It was found 
that non-linear Schrodinger equation aptly describes the 
pulse propagation in picosecond regime. Generalization 
of the NLSE is neccessitated to take into account higher 
order dispersion, self-steepening of the pulse due to the 
dependence of the slowly varying part of the non-linear 
polarization on time and the delayed effect of Raman 
response for describing optical pulse propagation in the 
femtosecond domain [4j. Higher order NLSE which cap- 
tures these effects is given by [!, Hj] , 



dz 

+ i(b 3 



2 dr 2 



6 dr 3 



3 ) = o. 
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It has both group velocity dispersion (GVD) and third 
order dispersion term, respectively represented by 02 and 
P$. Parameter b\ takes into account Kerr non-linearity, 
whereas self steepening effect arising from time varia- 
tion of the non-linearity (| <j> \ (f>) t has been incorporated 
through the coefficient 62. The terms with coefficient 63 
and 64 account for the self frequency shift, originating 
from delayed Raman response. 

This model, unlike NLSE, is not integrable in general. 
A few integrable cases have been identified; these are 
known as (i) Sasa-Satsuma case Q, (ii) Hirota case [1] 
and (iii) derivative NLSE of type I and type II A 
few restrictive special solutions of bright and dark type 
have been obtained for this system [lfj-[l2]]. Effects of 
these higher order terms on pulse propagation have been 
studied numerically quite extensively la. 1 1 31: some special 
solutions to this system are also known [14] . 



Third order dispersion becomes important for fem- 
tosecond pulses, when group velocity dispersion (GVD) 
is close to zero. It can be neglected for optical pulses, 
whose width is of the order of 100 femtoseconds or more, 
having power of the order of 1 Watt and GVD far away 
from zero. For short range propagation, 64 can also be 
ignored. The resulting equation can be written in a di- 
mensionless form by appropriately scaling variables: 



ip = n 



(f>, t 



and a; 



(2) 



where fi — [i x half width] of the incident pulse. The 
dimensionless equation obtained reads: 

.dip d 2 ip 2 , 
z- + ai — + H ^ 



+ i( a4 _(|^| 2 ^) + a5 -0 — (|^| 2 )) = 0, 



where. 



sgn(/3 2 



,o 4 



and 



a 5 
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In a previous study by some of the present authors, 
localized wave packets with nontrivial chirping were 
emerged as exact solutions of this equation |15|]. It was 
found that chirping has a kinematic component, deter- 
mined through initial conditions and a dynamical com- 
ponent, originating from self-steepening of pulse and its 
delayed Raman response. These two chirpings have in- 
verse characteristics and can be used for chirp control. 
Further, periodic solutions, with sinusoidal functions in 
fractional linear form, have also been identified. 

In this paper, we systematically investigate wide class 
of solutions of this system containing both bright and 
dark solitons. Interestingly, the equation governing dy- 
namics of intensity is found to be the one that governs 
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dynamics of field amplitude in non-linear Schrodinger 
equation with a source jl?], EH ■ This connection is then 
profitably used to obtain exact solitonic intensity waves 
of this higher order NLSE, possessing nontrivial chirp- 
ing behavior. The NLSE type bright solitons are found 
to propagate at superluminal velocity, with dark solitons 
complementing them in parameter domains, which are 
experimentally accessible. Interestingly, both the solu- 
tions also complement each other in velocity space. Fur- 
ther, we find that dark soliton which are static in systems 
like NLSE, are propagating in this system solely due to 
presence of higher order nonlinearities. These solutions 
show chirality, which depends on the sign of GVD param- 
eter. Depending on the kinematic chirp parameter two 
distinct domains are obtained. In the first case, if soli- 
tons are chiral for anomalous dispersion (a\ > 0), then 
they are non-chiral for normal dispersion (a± < 0). In the 
second case, solitons are strictly chiral in both the dis- 
persion domains. We then obtain more general localized 
solutions, through a conformal transformation, which can 
travel at both sub and superluminal velocities. In order 
to classify and gain physical understanding of the ob- 
tained solutions, we then use a pseudo-potential picture, 
which shows that the solutions separate into bounded 
periodic oscillations around the well minima and un- 
bounded motion, separated by the so called separatrix 
solution. The separatrix solution is identified and the 
connection of this system with hydrodynamical equations 
are pointed out. Interestingly, the localized kink type 
solitons of this system representing the separatrix can be 
mapped to the separatix solutions of the elliptic function 
equation, through the aforementioned conformal trans- 
formation. 

The paper is organized as follows. In the following 
section, we explore the localized dark and bright soliton 
solutions after connecting the present dynamical system 
with the NLSE. Chirality of the solutions and the effect 
of chirping have been systematically investigated. In Sec. 
Ill, we find more general solutions, through a fractional 
linear transform and show they can be both sub and su- 
per luminal depending upon various parameters. Section 
IV deals with conclusion and future directions of work. 



II. EXACT SOLUTIONS TO HIGHER ORDER 
NON-LINEAR SCHRODINGER EQUATION 
(HNLSE) 

We seek a travelling wave packet type solution of the 
HNLSE, which can generally be written as: 

ib(x,t)=p(0e lx{ Z ) - (5) 

Here the travelling coordinate £ = a(t — ux), and p and 
X are real functions of £, a is the scale parameter and 
u = 1/v, where v is group velocity of the ansatz solution. 
Substituting (0 in ©, and equating real and imaginary 



parts yields two coupled equations: 

/ I r> 2 II: 2 II 

~aup + La a\x p + a a\X P 

+ 3aci4p 2 p' + 2aa$p 2 p' — 0, (6) 
aux'p + a 2 a\p" — a 2 aix' 2 p 

+ a 2 p 3 - aa 4X 'p 3 = 0. (7) 
Equation ((6|) can be exactly integrated to yield: 

, u c (3q 4 + 2a 5 ) 2 

X - 7, 1 o a P ' W 

Zaai aaip z Aaa\ 

where c is a constant of integration, to be determined by 
initial conditions. Notice that the phase has a nontrivial 
form and has two intensity dependent chirping terms, 
apart from the kinematic first term, which is of usual 
e i{kx-wt) type. As is evident, the second term is also 
of kinematic origin and is infact present in the linear 
Schrodinger equation. The last term is due to higher 
order nonlinearities, which has a dynamical origin, and 
leads to chirping that is exactly inverse to that of the 
former. This is a novel form of self-phase modulation. 
Using the above expression, Eq. ([7]) can be written as: 

8i P " + 8 2 p + 8 3P 3 + 8 iP 5 = 4, (9) 

p A 

with 8 X = a 2 ai 2 , 8 2 = ( " 2+2a4 4 c+4Q5c) , 8 3 = (2ai ~" a4) and 

/) (a<i— 2a 5 )(3a4+2a 5 ) 

t>4 - is ■ 

Note that, in parameter regime, where 84 = 0, the am- 
plitude dynamics will be like the NLSE case, albeit with 
dressed parameters. Phase dynamics, however will not be 
the same as that of NLSE, and will depend on, whether 
(04 — 205) = or (304+205) = 0. In the latter case only, it 
will be exactly like NLSE, devoid of dynamical chirping. 
As shown in Ref. [15] . localized solutions of the above 
equation, in these parameter domains, show directional- 
ity and hence are chiral, and the propagation direction 
for these solitons is decided by sign of 04. Because, these 
solutions satisfy NLSE with dressed parameters, many of 
the parametric restrictions on solution space are relaxed, 
for example both dark and bright solitons are present in 
both anomalous and normal dispersion regimes. 

Defining k = 8 lP ' 2 + 8 2 p 2 + f p 4 + f p 6 + Eq. © 
can be written conveniently in terms of a (where a = 
p 2 = t/j*ip) to yield: 

d -la" + 28 2 a+ 3 la 2 + ^ = k, (10) 

where, k is to determined via initial conditions. Moti- 
vated by the similarity of the above intensity equation 
with that of amplitude equation of NLSE in the presence 
of a source [13, [3 1 below we consider different ansatz 
solutions. 

Case I: First we study localized soliton ansatz of the 
type 

<t = <to ^l + 7 sech(^)j , (11) 
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which allows for both dark and bright solitons, depending 
on the sign of 7. Consistency of the above solution yields 
the following conditions: 



6xa , 40 4 cr o s 



WW 
2 

2r] 2 



0, 



3 

A9 4 a Q 3 = 0, 

- 20 2 a o + 36W + 46W - 0, 



and 20 2 o"o H 1 fe = 0. 

2 3 

Solving first three of these relations, one finds 



30 3 2 IQO^i 
2 (90 3 2 



and 7 = - 



320. 



2<?4 ; 



(12) 
(13) 
(14) 
(15) 

(16) 
(17) 



Since, 0i > 0, reality of rj demands 04 > 0. Similarly, 
reality of 7 requires, 903 2 — 320204 > 0, whereas positivity 
of a leads to 03 < 0. These conditions, together with 
the above relations, lead to: 



(a4 — 205) > 0, u > — -, and 
(9a4 2 — 26) u 2 — (360104) u 



(18) 



+ (36ai 2 - 46c(a 4 + 2a 5 )) > 0; (19) 

where 6 = (04 — 2a 5 )(3a4 + 2a 5 ). 

In case of grey solitons, i.e., 7 < 0, for positivity of a 
requires 7 > — 1, which gives, 

(l5a4 2 — 26) u 2 ~ (600104) u 

+ (60ai 2 - 86c(a 4 + 2a 5 )) < 0. (20) 

We note that a\ can be positive or negative. Depending 
on the value of 04, u > shows possible chirality of 
the solutions. In case when 7 = — 1, intensity vanishes 
at origin and these solutions are known as dark solitons. 
Interestingly, derivative of phase (from Eq.®) becomes 
singular at that point, provided c ^ 0. This inturn im- 
plies that phase exhibits a discontinuous jumps across 
origin. Dark solitons of this kind have long been know 
in case of NLSE where they are strictly static [3j . Veloc- 
ity space of dark solitons in our case can be investigated 
using Eq. ([20| . One can easily check that these soli- 
tons in general are not static, unlike NLSE, and exhibit 
superluminal propagation with velocity depending on c. 
However, if (3a4 + 4a5)(a4 +4as) = 0, then these solitons 
can be static. Interestingly, in the limit 04,5 — > one 
can immediately see that Eg. (f2"0)) only allows for static 
solitons which corroborates with known results. This es- 
sentially implies that presence of higher order nonlinear- 
ities allows these kind of dark solitons to propagate with 
non-zero velocity. 



The last consistency relation (Eq.(|15p) determines u, 
the soliton velocity in terms of k: 



903 3 



2C3C4 



640 4 ''fc = 0; 

(M 2 +2a 4 c+4a 5 c) 



we note that 02 
contain u explicitly. As defined earlier, 

_/2 



and 03 



k = e x - 



03 n 04 r, C 



(21) 

(2ai — 1104) 
2 

, (22) 



4a 2 6 a 

is a constant of motion, and hence its value remains un- 
changed under time evolution. At £ — ¥ 00, it is seen that 
a' — > and a — > <tq, which allows one to write: 



k = 



3 2 3 27 3 3 8 4 c 2 
"8~0 4 _+ 512 07 ~ 3~03~' 



(23) 



Using this, along with Eq. (|2ip . yields an equation for u 
involving the chirp parameter: 



1350 3 4 - 360 2 3 2 6 + b 3 c 2 = 0. 



(24) 



When the constant of integration, c = one finds a 
quadratic equation for u: 



u 2 (15a 4 2 - 46) - 60aia 2 a 4 u + 60ai 2 = 0. 



(25) 



Reality of u can be ensured by demanding positivity of 
the discriminant in the above equation, which leads to 
04 > 0; this condition is trivially satisfied as it emerges 
from the reality of 77. 

We now systematically investigate the chirality of the 
solutions. For anomalous dispersion regime, a\ > 0, ve- 
locity is found to be positive and bounded from above, 
< v < t^-. For ai < 0, v > or v < -^-r. For 

2ai 1 5 — 2|ai| 

04 > 0, it is straightforward to see that both the roots of 
Eq. (1231) are negative, hence the possibility of v > is 
forbidden. Hence, the solution is unidirectional and show 
chirality, with the direction of propagation governed by 
the sign of <X\. Typical experimental values correspond- 
ing to pulse propagation around 100 femtosecond region 
are a\ — 0.5, a 4 = 10 _1 ' 5 ,a 5 = 10~ 2 . Using these values, 
one finds v = 0.01275, which after appropriate rescal- 
ing leads to v = 6.3 x 10 10 m/s, yielding superluminal 
velocity. We have checked that in the currently exper- 
imentally accessible parameter domain, all the solution 
are superluminal. 

Case II: Next we consider kink type solutions: 



a = a 1 



7tanh(— ) 
V 



(26) 



where | 7 | < 1 due to positivity of intensity. Proceeding 
like the earlier case, one finds four consistency conditions, 
three of which yield: 



(27) 
(28) 

(29) 





30 3 


^0 = 






80 4 ' 




160104 


v 2 = 


30 3 2 



and 7 = — 



1 (90 3 2 - 320 2 4 ) 
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Reality of 77 requires 84 < 0, which along with positivity 
of (To implies 9 3 > 0. This clearly implies that, the kink 
type solution and the earlier found soliton live in mutu- 
ally disjoint parameter and velocity spaces. Further, as 
in the earlier case, the reality of 7 requires, 

(9a 4 2 - 26) w 2 -(36aia 4 ) u+(36ai 2 - 46c(a 4 + 2a 5 )) > 0, 

(30) 

to be satisfied. Since (9a 4 2 — 26) > 0, the above inequal- 
ity holds for all real values of u, provided the discrimi- 
nant is negative. This leads to restrictions on the chirp 
parameter: 



c < 



18ai 2 



(a 4 + 2a 5 )(9a 4 2 - 26)' 



(31) 



When above relation is not obeyed, the allowed u values 
are forbidden between the two roots. As noted above, 
positivity of intensity demands I7I < 1, which translates 

I 



to: 



66» 3 2 - 326> 2 6» 4 < 0, 



imposing an upper bound on c: 

6 ai 2 



c < 



(a 4 + 2as)(3a 4 2 — 6) 



(32) 



(33) 



Since, -, — . „ 6 V-,„ — 2— rr > 0, for existence of these solu- 

' ((24 + 2(15 ) (3(I4 Z — 0) ' 

tions, c has to be a non-zero negative quantity. Like the 
earlier case, the consistency condition 0% > requires 



1 2ai 

v a 4 



(34) 



For the normal dispersion region, this condition gives, 
> v > ■ — TTi - ~n — and for anomalous region (ai > 0) 

2|ai|a,2 ° v 1 

above inequality leads to v < or v > 



<l4 



if 



(a 4 + 205) (3a 4 2 — 6) 



> c > 



6ai 2 



3a 4 3 + 2a 4 2 a5 — 12as 2 a 4 — 8a^ 



(35) 



When a% > 0, Eq. ([32)) forbids the negative values of 
velocity. Hence for this interval, we have strict chirality 
for the solutions. 

The remaining consistency condition yields an equa- 
tion for u: 



99 3 3 



2 6 3 6 4 - 646 4 2 k = 0, 



(36) 



which has the same form as in the previous case, albeit 
the value of k will be different. Constant of motion, k, is 
determined by evaluating it at £ = : 



3 02 



21_6^ 

128 6 A 2 



8 8 A c 2 



2 64 

Hence, equation for u reads 



3 e 3 



26o 2 



(37) 



1086> 3 4 - 726 2 6 2 b + 12(9 2 V + 6 3 c 2 = 0, (38) 

which is a quartic equation. We have not been able to 
find real roots of this equation for experimental values 
around 100 femtosecond pulse width. This can be proved 
by looking at the sub determinant structure of this quar- 
tic equation, i.e., Ai, A 2 , A3, and A 4 [H The 
sign of Ai (where i = 1,2, 3, 4), will determine the char- 
acteristics of the roots. Lets us take the list of signs of 
the Ai, i.e., [51,82,83,54], where Sj is the sign of Aj . 
The number of sign change in the sequence reveals the 
number of real pairs. If the sequence has n number of 
change in sign, then n pair of complex roots exist. Here 
Ai, A 4 > are always satisfied in the allowed domain 
i.e., 64 < 0. Hence for a real solution, we must have 
A 2 ,A3 > 0. In Fig. [T]we have plotted A 2 and A3 as 



a function of c and / respectively, one can see from the 
graph that both quantity do not share same sign for any 
value of c and /. 



III. CONFORMAL TRANSFORMATION 

Apart from ansatz solutions discussed above, below we 
show that Eq. fTUl) can be solved quite generally, using a 
conformal mapping involving Jacobi elliptic functions. 
Case I: When 64 ^ 

In this case, Eq. (|10[) can be written in a simpler form 
by going to y variable, y = a + 



■Ml,, 



y" +py + qy 3 



where p = ±(26 2 - f^), q = ^ and r = £( 

3 A 3 — k). In case when r = 0, above equation coincides 
with the equation satisfied by Jacobi elliptic functions. 
Below we show that, for some real constants A, B and 
C, the conformal transformation: 



r = 0, 



(39) 



()4«i" 



y = 



A + Bf 
l + Cf 



(40) 



solves equation ([39]) . with /(£, m) being one of the twelve 
Jacobi elliptic functions with modulus parameter m [13, 
li~8l ]. These elliptic functions satisfy: 



/" = af + bf 



(41) 



with appropriate real constants a and 6. For example, 
for / = cn(£, m), a — (2m — 1) and 6 = —2m. The claim 
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FIG. 1. A 2 ,3 are plotted above for ai = 0.5, a 4 = 0.03162, 
as = 0.01581 + / and c is chirping parameter 



hence can be solved exactly to yield: 

bA + aBC + aAC 2 
ABC 3 ' 
-3bA - 2aBC - aAC 2 



E a = 
V = 

q = 



C(-B + AC) 

-b - aC 2 



B(B - AC) ' 
bAB + bA 2 C + aB 2 C + aABC 2 
C 2 (B-AC) ' 



(46) 
(47) 

(48) 

(49) 



Therefore, this shows that, provided above relations hold, 
Eq. ([59"]) is indeed solved by (|3D]). Notice the depen- 
dence oi Eq on A, B and C, which simply shows that 
the initial values in Eq. (|4"Tj) do play a role in deter- 
mining A, B and C, and hence they can not be fixed, 
given p, q and r. Since, original equation parameters a* 
(i = 1, 5) can be expressed in terms of p, q and r, one 
can determine the parameter regime for a given solution 
using these relations. The solutions presented in Ref. 
[lH form a subclass of the ones found above. For m = 1, 
one gets localized solutions and are often of experimen- 
tal and technological interest. Setting m — 0, singles out 
periodic solutions, which were reported in Ref. [1 51]. 

For explicitness, we consider / = sech(£), which gives 
a = 1,6= —2 and Eq is found to be 0. The consistency 
equation becomes 



2A + BC + AC 2 = 0, 
C 



A{B - AC) 



q, and 



2AC + B 
~ B~ AC ~~ 

2A 2 
C(B - AC) 



p, (50) 

= r (51) 
(52) 



After eliminating B and C one gets, 



that expression (l4Tfl) solves Eq. (1591 can readily be seen, 

by defining the first integral Eq = 4j and 

substituting (gS) in Eq. Use of Eq. (gTJ leads to 

the following consistency conditions: 



r+pA + qA 3 =Q, (53) 

implying that velocity of soliton is governed by back- 
ground amplitude in this non-linear system. Since we 
know the values of p, q and r, by fixing the value of u, 
we can get the corresponding values of B and C in the 
form, 



- ABCE + 4AC 2 E + P A + qA 3 +r = 0, (42) 
aB - aAC + pB + 2pAC + 3qA 2 B + 3rC = 0, (43) 

- aBC + aAC 2 + 2pBC + pAC 2 + 3qAB 2 + 3rC 2 = 0, 

(44) 

bB - bAC + pBC 2 + qB 3 + rC 3 = 0. (45) 



Notice that these coupled algebraic equations are nonlin- 
ear in A, B and C, but are linear in p, q, r and Eq, and 



C 2 = ™Ll an d B 2 = ^ - 8 A 2 . (54) 
r qr 

For the experimental values, ai = 0.5, 04 = 0.0316228, 
05 = 0.01 and putting u — 10 5 (sub-luminal), we get 
A = 1.009 x 10 7 , B = - 6.51 x 10 7 and C = 6.7453. 
In the case of super-luminal propagation, u — 10, we 
get A = 1536.46, B = 1307.78 and C = 1.05128. As is 
evident, A, B oc u and C ~ 0(1). 

For / = tanh(£), in Eq. (gTJ| one finds a = -2 b = 2 
and Eq = 0.5. In this case, the consistency conditions 
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become 



-A + BC + AC 2 
BC 3 

-6A + ABC + 2AC 2 
C{-B + AC) 

-2 + 2C 2 



B(B - AC) 
-2AB 2 - 2A 2 C + 2B 2 C + 2ABC 2 
C 2 (B-AC) 



= 1, (55) 

= p, (56) 

= q, and (57) 

= r. (58) 



After simplification, we get a quadratic equation for u 
along with other consistency relations given by, 



9a 4 2 b 

u (— h — ) — 9aia4ii + (a,4 + 2a$)bc + bai = 0, 



90 3 3 - 4:89 2 9 3 9 4 - 646/ k = 0, 
46> 4 



C 2 = -^A 2 , 



B = 



A 
C' 



(59) 
(60) 

(61) 
(62) 



Once u is known, A can be evaluated from Eq. (|60[) 
through k dependence on A, then C and B can be ob- 
tained using Eq. (|6"Tj) and f[62]) respectively. For the 
consistency of Eq. (|6"Tj) , #4 < must be satisfied. Reg- 
ularity of solution requires C < 1. Also, the positivity 



of cr demands 



A- 1 



36> 3 



> 0, which gives restriction on 



l-C 

the velocity of the soliton. 

Above obtained solutions can be better understood us- 
ing following pseudo particle picture. Eq. (|39|) can be 
thought of as equation of motion for a classical particle 
of unit mass, with displacement given by y, moving un- 
der influence of a nonlinear force. The same also holds 
for Eq. (|4ip. albeit the nonlinear force does not have 
a constant term. Also notice that both these dynami- 
cal systems are of one degree of freedom and their phase 
spaces are two dimensional, only position y and momen- 
tum y' is required to describe the dynamics completely. 
Looked in this setting, mapping Eq. (1401) simply says 
how the two phase spaces are related; more precisely it 
shows how phase space for Eq. (l39t can be generated 
by knowledge of phase space for Eq. (|4"Tj) . Since, the 
mapping is conformal, the singularity structure of both 



the phase spaces are identical, modulo movable poles. In- 
fact, it indicates that Eq. (j39j) is integrable in the sense of 
Painleve, since Eq. (grj is integrable [H, The kink 
type solutions of the parent Eq. ([3]), in this setting, can 
be identified with the separatrix solutions of Eq. (|41[) via 
the conformal mapping 21]. Similarly, singular solutions 
of Eq. (|41|) can be seen to be related to unbounded solu- 
tions lying outside separatrix, and periodic solutions can 
be seen to be related to bounded solutions lying inside 
separatrix. 

It is very interesting to note that solutions to Eq. 
([3T)|) are actually solutions to modified KdV equation: 
v X xx + Piv 2 v x + p2V t — 0. This can be seen using the 
travelling variable £ = x — vt, integrating out the equa- 
tion once, and identifying p\ — Zq, p 2 = — f and r as 
the constant of integration. So, at the traveling variable 
level, which restricts one to only one soliton solution, one 
can say that the modified KdV solitons can be simulated 
by considering these higher nonlinearities in optical fiber. 
Case II: When 4 = 

In this case, the above analysis is not valid since the 
degree of the Eq. © changes and can be written as: 



p'a + q'a 2 



0. 



(63) 



with p' = ^4, q' = ^ an d r ' = -|^. The above equa- 
tion can be mapped to Eq. (|41|) via a transformation: 



(64) 



A + Bf 2 
C + Df 2 ' 



which is along the same lines as (|40|) . Following the same 
procedure as the former case, using the first integral Eq = 
4j 2i W e find the consistency conditions as: 

ABC 2 E - AACDEq + p'AC 2 + q'A 2 C + r'C 3 = 0, 

(65) 

AaBC 2 - AaACD - 12BCDE Q + 12AD 2 E + p'BC 2 
+ 2p'ACD + 2q'ABC + q'A 2 D + 3r'C 2 D = 0, (66) 

3bBC 2 - 3bACD - AaBCD + AaAD 2 + 2p'BCD + p'AD 2 
+ q'B 2 C + 2q ABD + ir'CD 2 = 0, (67) 
- b'BCD + b'AD 2 + p'BD 2 + q'B 2 D + r' D 3 = 0. 

(68) 

Again we observe that, the above equations are linear in 
p' , q' , r' and Eq, and hence can be solved to give: 
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_ bC 2 -2aCD 

~ — — ' (69) 

2(-3bBC 2 - SbACD + AaBCD + 2aAD 2 ) , , 

P = ~ D(AD - BC) (?0) 

6(-bC 2 + aCD) 

q AD-BC ' [ ' 

, _ -bB 2 C 2 - AbABCD + 2aB 2 CD - bA 2 D 2 + AaABD 2 

T ~ D 2 (BC-AD) ' ( ' 



This shows that, the Eq. (f63|) is an integrable equation, 
in the sense of Painleve, and the explicit solutions can 
be expressed in terms of Jacobi elliptic functions. By 
choice, of appropriate values of A, B, C and D one can 
find out the parameter regime in which the given solution 
is valid. Further, solutions to Eq. actually satisfy 

KdV equation: v xxx + p\W x + p2Vt = 0, in travelling 
variable £ = x — vt, with p' = —v/3, 2q' = a and r' being 
an integration constant. It is very interesting to note 
that the condition 64 — 0, can be fulfilled if 04 = 2a$ 
or 3a4 = — 2a§. In the case when the former is true 
then, the intensity profile of these solutions is exactly like 
NLSE solutions except with a non-trivial phase chirping. 
In the latter case, the solutions do not have this non- 
trivial chirping and the solutions are indistinguishable 
from NLSE solutions with appropriate coefficients. 



IV. CONCLUSION 

In conclusion, we have shown that nonlinear 
Schrodinger equation in the presence of self-steepening 
and self- frequency shift, possesses rich travelling wave dy- 
namics, with non-trivial chirping. Self-phase modulation 
revealed here indicates the possibility of effective con- 
trol of pulse dynamics through chirp management. The 
dynamics of the intensity is governed by the non-linear 
Schrodinger equation with a source. Hence, intensity 
waves show solitonic behavior, in complete analogy with 
the integrable NLSE. These solutions exhibit superlumi- 
nal propagation in experimentally accessible parameter 
domain. Further, these solutions reveal directional prop- 
agation, akin to hydrodynamical equations like KdV, 
with which the obtained solutions have close connection. 
We have found that dark solitons in this system can prop- 
agate with nonzero velocity, unlike in the case of NLSE. 
More general solutions obtained through fractional linear 
transform, showed both sub and superluminal behavior. 
We identified the separatrix solution bounding the reg- 
ular oscillatory motion from the unbounded ones. This 
strengthens the hope that this system may be integrable, 
since separatrix is not known to exist in a system with 
chaotic dynamics. 
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